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We have derived algorithms and techniques to precisely co-register laser altimeter profiles with gridded 
Digital Terrain Models (DTMs), typically derived from stereo images. The algorithm consists of an initial 
grid search followed by a least-squares matching and yields the translation parameters at sub-pixel level 
needed to align the DTM and the laser profiles in 3D space. This software tool was primarily developed 
and tested for co-registration of laser profiles from the Lunar Orbiter Laser Altimeter (LOLA) with DTMs 
derived from the Lunar Reconnaissance Orbiter (LRO) Narrow Angle Camera (NAC) stereo images. Data 
sets can be co-registered with positional accuracy between 0.13 m and several meters depending on the 
pixel resolution and amount of laser shots, where rough surfaces typically result in more accurate 
co-registrations. Residual heights of the data sets are as small as 0.18 m. The software can be used to 
identify instrument misalignment, orbit errors, pointing jitter, or problems associated with reference 
frames being used. Also, assessments of DTM effective resolutions can be obtained. From the correct 
position between the two data sets, comparisons of surface morphology and roughness can be made at 
laser footprint- or DTM pixel-level. The precise co-registration allows us to carry out joint analysis of the 
data sets and ultimately to derive merged high-quality data products. Examples of matching other 
planetary data sets, like LOLA with LRO Wide Angle Camera (WAC) DTMs or Mars Orbiter Laser Altimeter 
(MOLA) with stereo models from the High Resolution Stereo Camera (HRSC) as well as Mercury Laser 
Altimeter (MLA) with Mercury Dual Imaging System (MDIS) are shown to demonstrate the broad science 
applications of the software tool. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

Images and topographic data are fundamental for the explora- 
tion of planetary surfaces. Topography on planetary missions is 
typically derived by laser altimetry in the form of profiles acquired 
along the spacecraft ground track (Smith et al., 2000, 2011 ; Smith 
and Zuber, 2000). Alternatively, photogrammetric reduction of 
stereo images is used to produce gridded Digital Terrain Models 
(DTMs) (Preusker et al., 2012; Scholten et al., 2012; Oberst et al., 
2010b; Haase et al., 2012). The laser profiles benefit from highly 
accurate individual point measurements, while the DTMs benefit 
from a contiguous large-area coverage. While both data sets are 
typically in good qualitative agreement, detailed comparisons 
often reveal lateral and vertical offsets in the raw data sets, in 
particular, when approaching the resolution limits of DTM pixel or 
laser spot sizes. In this work we have refined the method to 
precisely co-register and determine offsets between the two data 
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sets (Glaser et al., 2010). These offsets offer valuable insights into 
various instrumental or operational effects, but may also be of 
great scientific interest. Also, the availability of precisely co- 
registered data sets offers great opportunities of joint analysis to 
improve accuracy, detail, and coverage of final data products. 

2. Previous work 

The Mars Global Surveyor (MGS) imaged the Martian surface 
with the Mars Orbiter Camera (MOC) and the Mars Orbiter Laser 
Altimeter (MOLA). Kim et al. (2000) used visible edges of surface 
features such as craters to co-register MOLA tracks with MOC and 
VIKING orbiter images. The proposed approach, however, can only 
find horizontal shifts between the two data sets if a crater can be 
detected in each of them. The accuracies of the shift also highly 
depend on the hierarchical and heuristic edge detection method 
that is used and are not mentioned in this work. Yoon and Shan 
(2005) showed a combined adjustment of MOC stereo images and 
MOLA data. Trajectory data, tie points of the images, MOLA ranges, 
and MOLA surface points are adjusted simultaneously in this work 




112 


P. Glaser et al. / Planetary and Space Science 89 (2013) 111-117 


using least-squares techniques. While an improvement for the 
ground location of MOLA surface points can be achieved, no 
condition is introduced in the algorithm minimizing the displace- 
ment between the stereo DTM and the MOLA profiles. Kolb and 
Okubo (2009) introduced a program to help co-registering MOLA 
data with Mars images obtained by MOC, the High Resolution 
Imaging Science Experiment (HiRISE) and the Context Camera 
(CTX). The user has to apply lateral shifts between images and 
MOLA tracks, but cannot align the two data sets in the 3-D space. 
No information on the fit accuracies can be retrieved in this 
approach. Di et al. (2012) and Lin et al. (2010) amongst others 
use surface matching techniques to register laser derived DTMs 
with stereo DTMs. Di et al. (2012) show surface matching using 
Chang'E-1 stereo images and crossover corrected laser DTMs (LAM) 
whereas Lin et al. (2010) match MOLA DTMs with High Resolution 
Stereo Camera (HRSC) and HiRISE DTMs. In both approaches the 
laser DTMs will contain tracks suffering from positional errors and 
interpolated data gaps. 


3. Data sets 

Several data sets were used to demonstrate the broad science 
applications of the software tool. 

3.1. LOLA 

The Lunar Reconnaissance Orbiter (LRO), launched in 2009, is 
equipped with the laser altimeter LOLA (Lunar Orbiting Laser 
Altimeter) and the camera system LROC (LRO Camera) featuring a 
pair of narrow- and one wide-angle camera (briefly termed NAC - 
Narrow Angle Camera and WAC - Wide Angle Camera). LOLA 
generates single laser pulses transmitted through a diffractive optical 
element, which splits up the laser pulse into five individual beams. 
For each beam the time of flight (range), pulse spreading (surface 
roughness) and the ratio of transmitted/returned energy (surface 
reflectance) are measured. From the pattern of the five laser spots, 
surface slopes not only along but also across the orbiter's track can be 
determined. From LRO's nominal orbit altitude of ~ 50 km, LOLA 
covers a ~ 50 m wide swath and the LOLA laser spot sizes are 5 m in 
diameter. At the laser pulse rate of 28 Hz, along- and across-track 
spacing of the spots is 10-12 m. The nominal ranging precision of the 
instrument is 10 cm and its vertical precision is < 1 m (Smith et al., 
2011). In this work we used parts of a LOLA track from orbit 3508 
(DOY 090 2010, primary mission phase). 

3.2. LROC NAC 

The two NACs, NAC-L and NAC-R, as part of the LROC system, 
are line scanners (5000 pix/line), operated in push-broom mode. 
From the nominal orbit height of 50 km, the two NAC cameras 
cover a 5 km swath with a length of ~25 km at an image scale of 
0.5 m/pix. For the production of DTMs, multiple observations from 
different orbits and appropriate convergence angles are required. 
Typically, one nadir-looking image is combined with an off-nadir 
image taken from the adjacent orbit with the spacecraft being 
tilted up to 20°. 

For our analysis, two NAC-L images (0.5 m/pix image resolu- 
tion) from subsequent orbits 1943 and 1944 (DOY 322 2009, 
primary mission) of the Apollo 14 landing site are used. The image 
pair NAC-L/NAC-L overlaps by 10.3% under perspective offsets, 
where stereo processing is possible. The solar incidence angle is 
58° for the image from orbit 1943 and 57° for the image from orbit 
1944 (At % 2 h). The DTM is computed at a grid size of 2 m, which 
is about half the LOLA spot size. The processing of the DTM was 
carried out by an adapted version of the DLR photogrammetric 


processing system, which was originally implemented for data 
from the Mars Express High Resolution Stereo Camera (HRSC) 
(Gwinner et al., 2009, 2010). 

3.3. LROC WAC 

The Wide Angle Camera (WAC) on-board LRO provided near- 
complete coverage of the lunar surface in stereo. The resulting 
DTM, GLD100 (Scholten et al., 2012), has a grid size of 100 m and 
covers the Moon from -79° to +79° latitude. 

3.4. MESSENGER MDIS and MLA 

The MESSENGER (MErcury Surface, Space ENvironment, GEo- 
chemistry and Ranging) spacecraft has completed three Mercury 
flybys and entered Mercury orbit in March 2011 (Solomon et al., 
2011 a, b). We have analyzed tracks of the Mercury Laser Altimeter 
(MLA) from the first Mercury flyby (Zuber et al., 2008) and 
matched them with a DTM derived from stereo data of the 
Mercury Dual Imaging System (MDIS) obtained during the second 
and third flybys of the spacecraft (Preusker et al., 2010). The DTM 
pixel resolution is 1 km in contrast to the MLA spot size which 
varies between 23 and 134 m during the flyby (Zuber et al., 2008). 

3.5. MGS MOLA and MEX HRSC 

The Mars Global Surveyor (MGS) entered Mars polar orbit in 1997. 
The on-board Mars Orbiter Laser Altimeter (MOLA) was in operation 
until June 30, 2001 (Neumann et al., 2002) and has obtained 600 
million range measurements (Neumann et al., 2003) with a spot size 
of ~75 m (Neumann et al., 2003). Mars Express (MEX) was launched 
in June 2003 and began its commissioning phase in January 2004 
(Chicarro, 2005). The on-board High Resolution Stereo Camera 
(HRSC) is a nine-look along- track stereo camera system. The DTM 
used in this work has a resolution of 100 m per pixel. 


4. Method 

The registration of the two data sets is performed by minimiz- 
ing the height differences of the n laser shots and the underlying 
DTM 

n 

^ ] (hdtrrii ~ hlaseri) — »min (1) 

i = 0 

with h laser . being the i-th height value measured by the laser and 
h dtm . the corresponding height in the DTM. A grid search based 
algorithm shifts the whole track from its nominal position 
in a user-defined search window, e.g. / = + 50 pixels in line and 
s = ± 50 pixels in sample direction and extracts the height value at 
each position from the DTM using a bilinear interpolation. 
At each shifted position the standard deviation is determined by 

2/*= o fidmiijs ~ h laser i) ^ 

n^l ( j 

The best match between the data sets at pixel size is found at the 
minimum of Eq. (2). By dividing the pixel position in a user- 
defined sub-pixel grid, e.g. step size of 1/30 pixel in line and 
sample direction, and using Eq. (2) again, a sub-pixel position of 
the best match can be found 

_ /ZF= o (hdtm u+subl ' S+subs ~ hlaser J 2 

^^/ + sub ( ,s + subs y Yl 1 ^ 1 

From the sub-pixel position, translation parameters x,y,z in 
sample, line, h direction or longitude, latitude, h respectively can 
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be extracted. In a second step, we wish to determine the 
accuracies of these parameters by means of least-squares techni- 
ques using the Gauss-Helmert model. The functional model is 

F(L,Xq) = Li — /(x+ Axo,y + Ay 0 ) — (z+ Azo) = 0 (4) 

Li is the observation vector containing the laser heights, /(x+Ax 0 , 
y + Ay 0 ) is a bilinear interpolation of DTM heights at the laser 
position and z+Az 0 is the vertical offset between the two data 
sets. To the translation parameters (x,y,z) that were derived by the 
grid search approach, we add the initial values Ax 0 ,Ay 0 ,Az 0 
representing the unknown vector X 0 in the least-squares approach. 
The initial values are all set to 0. The condition is that the difference 
of each laser height and the interpolated DTM height at the laser spot 
position corrected for the offset should be 0. 

In order to give accuracies for the resulting fit of the laser track 
with the DTM, the variance covariance matrix of the adjusted 
unknowns, X xx , has to be evaluated: 

^xx — SoQxx (5) 

On the main diagonal we find the squares of <r Axo , o - Ayo » ° Az 0 pointing 
in sample, line and h direction (longitude, latitude, h ) respectively. 
If the result of the grid search is correct then the resulting Ax 0 , Ay 0 , 
Az 0 should not be significant, therefore numerically smaller than 
their corresponding accuracy levels c A * 0 , o- A jv °Az 0 - 


5. Results 

LRO LOLA-NAC : A sample analysis for co-registering a LOLA 
track (243 laser spots) with an Apollo 14 landing site DTM is 
shown in the following. The nominal position of the LOLA track 
from orbit 3508 and the corresponding DTM heights are plotted in 
Fig. 1. Significant differences in profile heights are noted (standard 
deviation of residuals: approximately 2.6 m). However, the differ- 
ences are systematic (probably due to pointing or orbit errors) and 
can be corrected. Small horizontal offsets of a laser track from its 
correctly matched position can result in substantial vertical offsets 
to the reference topographic model when a laser track passes 
through rough terrain, while vertical differences on flat terrain are 
comparably small (see Fig. 1). Craters constrain the horizontal and 
vertical shifts and are ideal for co-registration. The first step of the 
grid search algorithm returns an integer pixel shift, see Fig. 2, left 
plot. Around this integer pixel position a sub-pixel grid search 
with a step size of 1/30 pixel is performed, see Fig. 2 right plot, for 
refinement. Applying the translations and re-plotting the two 
profiles, the standard deviation of the height differences drops to 
21.22 cm, see Fig. 3. The height differences vary within an interval 
of + 50 cm, which probably originate mostly from differences in 
the effective height resolution of the DTM and the laser profile. 

Least-squares result : The standard deviations of the adjusted 
unknowns are <t Axo = 8.75 cm, o- Ayo = 8.71 cm and o- Azo = 1.60 cm 
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Fig. 1. The nominal LOLA profile and the corresponding DTM heights. At the top of 
the plot the difference between the two profiles is shown. Height values are given 
in meters with respect to a sphere with R= 1737.4 km. (For interpretation of the 
references to color in this figure caption, the reader is referred to the web version of 
this article.) 


Latitude 

Fig. 3. DTM heights and the corresponding LOLA profile after applying the final 
translation parameters. The misplacement between the two data sets is minimal. 
(For interpretation of the references to color in this figure caption, the reader is 
referred to the web version of this article.) 
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Fig. 2. Left: residual field of integer shifts of the grid search algorithm. The minimum is well defined. Right: the result of sub-pixel grid search. (For interpretation of the 
references to color in this figure caption, the reader is referred to the web version of this article.) 
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(see Section 4), which attest to the robustness of the fit. The 
configuration of the NAC DTM with 2 m resolution and the LOLA 
laser track with its 5 m spot size results in a sub-pixel accuracy 
after the co-registration for both data sets. 

LRO LOLA-WAC : The shift needed to co-register the LOLA profile 
from orbit 10,601 with a subset of the GLD100 is adequate for the 
entire profile. Even though the track extends over 300 km, there 
are no obvious rotational offsets, scale differences or trends. The 
standard deviation of the height differences after the registration 
is 19.12 m, which is less than 1/5 DTM pixel. 

The main problem when matching DTMs with relatively large 
DTM pixels and laser profiles with high spot to spot densities is that 
multiple, accurately measured laser height values are compared with 
similar or even the same interpolated, smoothed DTM pixels. In the 
example shown in Fig. 4 it can be seen that ~ 7 LOLA spots fall within 
the same DTM pixel and therefore will be compared to almost the 
same height value retrieved by the interpolation in the DTM, which 
introduces high correlation in the least-squares adjustment. For 
this purpose we co-registered 3 versions of LOLA profile 10,601 to 
the GLD100 resulting in a standard deviation of height residuals of 
19-20 m. First we took the original profile with all 25,920 spots, 
which lead to g Ax = 2.47 m, G Ay = 1.54 m, g Az = 6.11 m. For the 
same profile three spots have then been averaged resulting in 8642 
spots and g Ax = 5.33 m, o Ay = 0.71 m, g Az = 2.89 m. Finally aver- 
aging the 5-spot pattern of LOLA to 1 spot and skipping every other 



Fig. 4. Distribution of LOLA points on the 100 m-WAC-DTM GLD100. About 7 LOLA 
spots (red dots) fall into the same DTM pixel (black squares) which introduces high 
correlations between the height observations of the laser and the interpolated DTM 
heights. (For interpretation of the references to color in this figure caption, the 
reader is referred to the web version of this article.) 


lead to 2611 spots and g Ax = 3.23 m, a Ay = 4.77 m, g Az = 0.38 m. 
The horizontal accuracy of the fit is best for the original profile 
(o Axy = 2.91 m) and is degrading the more the profile is averaged 
{o Axy = 5.37 m, a Axy = 5.76 m). On the other hand the high spot 
density of the original profile in comparison to a GLD100 pixel leads 
to inferior vertical accuracy of the fit compared to the averaged 
profiles, which is described by the correlation in the least-squares 
adjustment (Fig. 5). 


Table 1 

Overview of comparisons shown in this work. The first column lists the pixel-sizes 
of the DTMs (D) used in this work, the second column shows the size of the 
footprints of the laser altimeters (I). The third column is the ratio of laser spot size 
and the DTM pixel. The forth and fifth column show the standard deviation of 
height differences between the two data sets after co-registration and the standard 
deviation of the fit. In the last two columns we show the ratio with the standard 
deviation of the fit. 


Missions 

D (m) 

L(m) 

L/D 

(-) 

h 

(m) 

dfit 

(m) 

L/dfit 

(-) 

D/a fit 
(-) 

LOLA-NAC 

0.5 

5 a 

10 

0.21 

0.16 

31.25 

3.12 

LOLA-WAC 

100 

5 a 

0.05 

19.07 

6.76 

0.74 

14.79 

MLA-MDIS 

250 

28-35 b 

0.112-0.14 

79.18 

50.37 

0.55-0.69 

4.96 

MOLA-HRSC 

100 

75 c 

0.75 

50 

13.43 

5.58 

7.45 


a Smith et al. (2010). 

b The distance from the MESSENGER spacecraft to the surface of Mercury for 
the shown example (see Fig. 7) varied from 440.050 km (36.005°N) to 346.979 km 
(53.9905°N). With a beam divergence of 80 prad for MLA, the resulting laser spot 
sizes therefore vary between 35 m and 28 m. 
c Neumann et al. (2003). 
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Fig. 6. MLA track matched with an MDIS DTM, computed from stereo images 
obtained during different flybys. (For interpretation of the references to color in 
this figure caption, the reader is referred to the web version of this article.) 




Fig. 5. LOLA track 10,601 on the GLD100. The height differences of the 300 km long 
profile, containing 25,920 laser shots show excellent agreement without obvious 
residual offsets and trends. (For interpretation of the references to color in this 
figure caption, the reader is referred to the web version of this article.) 


Fig. 7. MLA laser track matched with an MDIS DTM from the MESSENGER orbital 
mission. The good match of the data sets suggests that altimeter, image, and navigation 
data are consistent. (For interpretation of the references to color in this figure caption, 
the reader is referred to the web version of this article.) 
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Fig. 8. Shifted MOLA track #12811 and HRSC DTM. (For interpretation of the 
references to color in this figure caption, the reader is referred to the web version of 
this article.) 



Fig. 9. Shifted LOLA track 1619 and NAC DTM (Orbits 598, 599). Jitter in the DTM 
is clearly visible in the resulting height differences and a frequency of ~ 1.7 Hz can 
be observed. (For interpretation of the references to color in this figure caption, the 
reader is referred to the web version of this article.) 


6. Other missions 

For other planetary data sets the relationship between DTM 
grids and laser profiles is usually quite different than for LRO data 
regarding laser shot spacing, laser spot sizes, DTM pixel sizes, and 
data accuracy (see Table 1). 

MESSENGER MLA-MDIS: An MLA track from the first Mercury 
flyby (Zuber et al., 2008) is matched with an MDIS DTM 
obtained during the second and third flybys of the spacecraft 
(Preusker et al., 2010) under difficult conditions, such as dras- 
tically changing spacecraft range and off-nadir pointing (see 
Fig. 6). With a laser spot size of ~1/10 of a DTM pixel the 
geometric relationship between the two data sets is similar to 
the LOLA-WAC example described earlier. The MLA track from 
the first flyby consists of 2360 laser points. After matching with 
the DTM, the standard deviation of the height differences is 
311 m in contrast to 349 m at the nominal position. A slope of 
~0.04° between the two profiles is identified, suggesting that 
probably errors existed in the navigation data of the early flybys. 
More recent DTMs and MLA tracks from the orbital mission phase, 
processed with more accurate knowledge of orbital position and 
attitude, lead to smaller residual offsets. The displayed MLA track 
(see Fig. 7) was recorded in April 2011 and cross-cuts crater 
Rustaveli. The accuracy of the fit is o Ax = 33.88 m, o Ay = 37.19 m 
and o Az = 2.43 m with a standard deviation of height residuals of 
79.18 m. 

MGS MOLA-MEX HRSC : Another example is given by matching 
MOLA tracks with HRSC DTMs. Our sample HRSC DTM, covering 
parts of the area of Argyre Planitia below 40° south, is very 
rough and has a high dynamic height range. The example (Fig. 8) 
shows parts of the shifted MOLA track from orbit 12,811 on the 
DTM, with which 1387 shots were co-registered. Given a laser 
spot size of ~75 m (Neumann et al., 2003) and a DTM pixel size 
of 100 m, the geometric relationship between DTM grid and 
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Fig. 10. Left: co-registered positions of LOLA spots plotted onto a NAC orthoimage. Right: laser spots color-coded for surface roughness, derived from pulse width. (For 
interpretation of the references to color in this figure caption, the reader is referred to the web version of this article.) 



Fig. 11. Left: Saddle Rock identified in the Apollo 14 NAC image near Cone Crater. Red square shows the area presented in Fig. 10. Right: Panorama of Saddle Rock (Hasselblad 
photographs AS14-68-9449/50/51) taken by Edgar Mitchell. The rock was described as “Volkswagen-size”. ©NASA/GSFC/Arizona State University/astroarts.org. (For 
interpretation of the references to color in this figure caption, the reader is referred to the web version of this article.) 
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laser spots is different to the examples shown before. The 
standard deviation of the height differences improved from 
65 m before to 50 m after matching, which is half the pixel size 
of the DTM. The accuracy of the fit is <r Axo =9.1 m, o- Ayo = 9.78 m 
and a Azo = 1.35 m. 

Review of the results: In Table 1 we summarize the results of 
this study. It is comprehensible that if the ratios of the laser 
footprint and DTM pixel size to the accuracy of the co-registration 
( L/ctfit and D/o^ f ) shown in the last two columns are >1, accurate 
fits can be made (at sub laser footprint and subpixel level 
respectively). The number of measurements and the accuracy of 
the data itself is very crucial but as a rule of thumb we can derive 
that larger ratios of L/D result in accuracies in the sub pixel level 
and laser footprint scale. We find that laser spots greater than 75% 
(see MOLA-HRSC and LOLA-NAC) of the DTM pixel are favorable 
for co-registering. 



Fig. 12. Offsets of 43 crossover improved LOLA tracks to an Apollo 17 NAC DTM 
with 1.5 m pixel resolution. The error bars (magnified by a factor of 10) vary from 
15 cm to 4 m with an average of ~ 50 cm. The majority of the offsets ( ~90%) scatter 
around the mean in a circle of ~ 15 m radius, which is the total position accuracy of 
the crossover improved LRO orbits (Mazarico et al., 2012). (For interpretation of the 
references to color in this figure caption, the reader is referred to the web version of 
this article.) 


7. Applications and future work 

A correct relative position of the laser track on the DTM is needed 
for many different applications in planetary mapping and science. 

• A comparison of the elevation profiles retrieved from a DTM 
and a laser profile may be used to assess the quality of the 
individual data products. Effects like image jitter (see Fig. 9), 
erroneous DTM trends (see Fig. 6), perhaps due to errors in the 
relevant navigation data, can be detected and corrected for. 

• After co-registration with DTMs, direct comparisons can be 
made between LOLA and orthoimages (which are typically 
derived from the same DTMs) at the limits of data resolution. 
In particular, with the correct laser spot position, roughness 
derived from the laser pulse width (Neumann et al., 2003; 
Smith et al., 2001; Garvin et al., 1999) can be compared with 
the orthoimage, and a calibration of laser pulse width for 
correct surface roughness can be made. The derived roughness 
can also be compared to the actual surface morphology in the 
image in full detail (see Fig. 10). Note the correlation of color 
coded laser pulse spread with surface morphology, the two red 
data points (5 m diameter) indicate high roughness as seen by 
the laser. The lower red colored laser point marks a rock; most 
likely this is ‘Saddle Rock’ visited by the Apollo 14 astronauts 
(Alan Shepard and Edgar Mitchell) (see Fig. 11). 

• Offsets between a camera derived DTM and laser tracks of the 
same orbiter can be indicative for misalignment between the 
two instruments. We detect persistent offsets between LOLA 
tracks and the NAC DTM, which are a strong indicator for small 
residual errors in the nominal pointing kernels (see Fig. 12). 

• A comparison of the topographic profiles of a DTM and a laser 
can be used to assess the effective height resolution of the DTM. 
While the DTM grid size alone does not necessarily represent its 
own effective resolution, a comparison of the DTM with the laser 
track, considered highly accurate, is helpful. By introducing a 
priori standard deviations for the laser and DTM heights, 
statistical tests (x 2 ) can be used to test if chosen parameters 
are appropriate. Hence, a statistically verified standard deviation 
of DTM heights can be given (Oberst et al., 2010a). 

• Co-registration offsets may be used to identify map-tie pro- 
blems, i.e. inconsistencies in the coordinate system parameters 
or erroneous assumptions in planet rotational model. A first 
analysis of MDIS DTMs suggests that comparisons with multiple 
laser tracks taken at different times can provide an observation 
for the librational motion of Mercury. Periodic shifts and offsets 
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Fig. 13. Connecting ridge - Candidate ESA Lunar Lander landing site. Left: gridded LOLA DTM at nominal position of tracks. Note that individual tracks are obviously slightly 
offset. Right: gridded LOLA DTM at co-registered position. (For interpretation of the references to color in this figure caption, the reader is referred to the web version of this 
article.) 
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have been detected by the introduced algorithm (Stark et al, 
2012 ). 

• Correctly co-registered laser tracks may be introduced to help 
improve terrain models, in particular to redo the bundle 
adjustment of the images. The registered laser track spots will 
serve as control points at fixed lateral and vertical coordinates. 
Likewise, misplaced laser tracks can be identified, corrected, or 
removed to yield improved gridded laser terrain models. Such 
work is currently in progress within our team to generate lunar 
landing site DTMs for the European Space Agency (ESA) (see 
Fig. 13; Carpenter et al., 2012; De Rosa et al., 2012). 


8. Summary 

It could be shown that the algorithm finds a best fit with sub- 
pixel accuracy between laser profiles and DTMs. Further a reason- 
able amount of laser shots is required to have a statistically reliable 
result. The use of this routine is not limited to one specific data set 
like LRO but is kept fairly independent of required data formats 
and missions. The range of application is big and the list we 
provided is most certainly incomplete. An overview of DTM pixel 
sizes, laser spot sizes, and the resulting standard deviations of 
height differences after matching is given in Table 1. 
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